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ABSTRACT 



There is a need for a fast method of analyzing bathythermograph 
traces and this need is approached by the means of a high-speed digital 
computer. The theory behind the computer program is outlined. Both 
synthetic and real data cases are run as examples using both data card 
decks and magnetic tape inputs. The program has been designed to read 
digitized bathythermograph traces and then analyze them objectively by 
Gaussian and non-Gaussian methods for the top, center, and base of the 
main thermocline. Additionally, such features as multiple thermoclines , 
inversions, and thermal transients are identified also and their key 
points are included in the informational data printout. 
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1 . 



Introduction. 



Since the introduction of the bathythermograph by Rossby and 
Montgomery (1934), and its further modification and perfection by 
Spilhaus (1938), there has been almost continuous thermal structure 
information funneled into the various collection agencies which has 
been handled manually under procedures outlined by LaFond (1951). 

This data collection has resulted in vast numbers of bathythermograph 
slides that have been noted, and then stored away for future analysis. 
In more recent years bathythermograph data has been stored on aperture 
cards, using techniques similar to those suggested by Sauer (1964). 

The Canadian Oceanographic Data Center, Scripps Institution of Oceano- 
graphy, and the Fleet Numerical Weather Central, Monterey, California, 
have further modified this procedure by a method similar to that sug- 
gested by Sauer and Hope (1967), where a device has been invented which 
will take incoming bathythermograph slides and convert them directly to 
digitized data format. This system described by Sauer and Hope (1967) 
consists of taking the grid/slide combination and transferring it to an 
aperture card by means of the BT Processor, a modified 3M Processor 
Camera. It has the advantage of preserving the original analog trace. 
The aperture card is then fed through a hopper device directly into a 
Caps Jeffree projector with a D-MAC curve-follower unit. The time re- 
quired for this unit to process one bathy thermogram trace is approxi- 
mately 45 seconds with a claimed accuracy of 1.0% of the full depth 
range and 0.1 degrees Farenheit. After exhaustive testing at the 
Naval Hydrographic Office as described by Stewart (1963) Sauer and 
Hope (1967) of the Canadian Oceanographic Data Center felt this degree 
of reliability in fact was not approachable. Accordingly, they 
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converted their system to record temperature to a minimum of every 0.1 de- 
grees Celcius and to the nearest whole meter. 

Papers by Tully (1961, 1964) show that there are certain features 
which must be identified in order to classify a bathythermograph trace. 
Tully states that these few parameters can be used to identify the thermal 
structure in any area of the ocean. Among these features is a potentially 
isothermal layer, which may or may not contain transients. There is a 
seasonal thermocline, in which the gradients may be regular or irregu- 
lar. Below this layer is a sub-thermocline layer in which the gradients 
are less than in the thermocline. Superimposed on these features are 
various types of short term fluctuations which possess varying degrees 
of significance, such as temperature inversions. 

Tully (1964) suggests in this same paper that the thermal structure 
can be defined by the following features: 

(1) Surface Temperature . 

(2) Potential Layer Depth . This figure includes the average 
depth to the top of the thermocline, plus any modifications due to inter- 
nal waves, etc. 

(3) Depth to the Bottom of the Thermocline . 

In many oceanic areas this depth has an associated distinct tempera- 
ture minimum and therefore is very definitive. In areas like the 
tropics it is usually not a significant feature. 

(4) Thermal Magnitude of Thermocline . To researchers this 
may be a very significant feature, because it can signify a boundary 
to restrict fish movement, water mass mixing, etc. 

(5) Other Less Definable Features . There are other parameters 
that may be considered, like the location of temperature minimum in a 
thermal inversion, mixed layer depth, etc. These features cannot be 
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fitted easily into a generalized machine processing format for ready 
identification. They are not common features, and may have special 
variables that depend solely upon the area of the ocean being investi- 
gated, or upon the season of the year. 

As Boston (1966) states, "until the information on bathythermograph 
slides can be conveniently analyzed by machine, much information about 
the thermal structure of the upper layers of the ocean is quite likely 
to remain hidden in an ever growing mountain of smoked and/or gold film 
glass slides, photographs of same, or punched cards." Two distinct 
problems arise: (1) the digitizing of the bathythermograph data to a 

form that can be processed by high-speed computers and data handling 
systems; and (2) specification of a program that will take this digi- 
tized data and process it for use in oceanographic research; the pro- 
gram needs the generality to be able to handle any shape trace, and flexi- 
bility so that further routines may be added or subtracted at will with- 
out destroying the basic integrity of the program. As we have seen, 
some progress is being made in the handling of the first problem. 

The second problem is the topic of this thesis; the objective 
analysis of digitized ba thy thermogram traces. The standard features 
as suggested by Tully (1964) are identified and printed out in a concise 
informative way. One or two additional features, like the identifica- 
tion of thermal inversions and transients, are included in an attempt 
to show that further features may be added to this program according to 
the requirements of the individual user. 
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2 . 



Method of Analysis of Thermal Structure , 



There will be two separate analytical methods developed for handling 
the objective analysis of the thermal structure. In his technical paper, 
Boston (1966) suggested the following two methods, on which this thesis 
is based. One method consists of a Gaussian temperature structure being 
assumed for the thermal parameter discussed. The second method is a 
common finite difference method, utilized to find inflection points in 
the various trace slopes. 

A prerequisite to both methods is that the data from the original 
analog bathythermograph trace be reduced to an ordered set of depths 
read at some predetermined temperature interval. From this ordered set, 
the thermal structure can be defined objectively in terms of certain pri- 
mary and secondary features. Both methods require only the first and 
second finite differences with depth; therefore only data from depths 
corresponding to some constant thermal interval (conveniently 0.5 to 1.0 
degrees) are required to handle the analytical technique. However, it is 
convenient to include temperature data corresponding to these points, as 
in the computer output there are locations where information as to the 
thermal magnitudes are as important as the spacial magnitudes. 

Thermal structure, as defined in this thesis, is the relationship of 
the temperature, T, to depth, Z, as a function of time and position. 
Furthering earlier ideas, Tully (1964), Tully and Giovando (1963), and 
Uda (1963) show in their papers that the thermal structure changes with 
season, latitude, and even oceanic region. However, Tully (1964) also 
suggested that the number of different features for a given application 
is small. Section A considers those features common to most oceanic 
regions within the upper 300-400 meters of the ocean. Section B considers 
those features that are of major concern in such fields as antisubmarine 
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warfare, biological forecasting, and fisheries that are not discussed 
or included in Section A. 

A. Primary Features: 

The primary features as suggested by Tully (1964) and Boston (1966) 
consist of several convenient parameters; (1) surface temperature, (2) 
depth of potential layer depth (upper isothermal), (3) depth of the 
bottom of the thermocline, and (4) magnitude of the thermocline. 

The potential mixed layer is that as defined by H. 0. Pub SP-105, 
Vol 5, and is considered to be isothermal, although this "isothermal 11 
region accepts small positive or negative gradients up to 0.3 degrees 
Farenheit per hundred feet. The termination of this nearly-isothermal 
layer is usually abrupt enough to clearly mark the region where the 
thermocline or other feature takes place. The depth of the mixed layer 
will be defined as the top of the thermocline. Likewise, the 

base of the thermocline, can be considered the top of the lower nearly- 
isothermal area. 

Once the top and bottom of the thermocline have been defined, with 
"TtL as the temperature at 2^ , and Tj-^as the temperature at 2 bt * 
the definition of the magnitude of the thermocline is 

AT = Tt -Tw, (1 > 

B. Secondary Features 

The secondary features are usually considered as time varitions in 
the primary features of Section A, and can occur at any point in the 

vertical water column being considered. 

1. Transients 

The term transient will be defined here as referring to short term 
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fluctuations of the thermal structure. It is usually the result of wind 
mixing, and vanishes over a period of a day or more. The program pro- 
duced for this thesis will recognize these transients, assuming that 
these thermal transients take the usual form of small but sharp gradients 
that appear as irregularities on the main thermal structure. Transients 
are often found in association with the well known "afternoon effect.” 

This program will aid in the verification of the prediction methods out- 
lined by Tabata and Giovando (1962). 

2. Multiple Thermoclines 

Heating and wind mixing conditions can also produce more than one 
stable thermocline. The growth of two stable thermoclines has been well 
documented by Tully and Giovando (1963). The occurrence of these multiple 
thermoclines is sufficiently common and vital to underwater acoustic pre- 
diction systems that they merit identification. 

3. Inversions 

Inversions are a common oceanic feature, an example of which has 
been described by Roden (1964), and must be identified if one is to ac- 
curately predict underwater sound transmission. The thermal inversion 
takes form when the temperature below the main thermocline (below Zfc>t) 
decreases with depth to a minimum, then increases markedly, and finally 
decreases again, but now usually at a more reduced rate than the first 
decrease. These thermal features are of extreme interest in underwater 
communications as the point of minimum temperature in this thermal 
structure forms the axis of a "sound channel” formed by the thermal bend- 
ing of acoustic rays by this inversion. 

The following methods will be used to identify these primary and 
secondary features within the computer program designed for this thesis. 
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c. 



The Gaussian Thermocline 



A gaussian, or normal, distribution (see Figure 1) of temperature, 
T, as a function of depth, Z, is given by 



A change in the numerical value of m causes a displacement of the curve 
in the vertical direction, Z, but does not alter the form of the curve. 

A change in CT has the effect of altering the scale of both coordinate 
axes. The smaller the value of CT , the more concentrated the curve is 
about the point Z = m. If the water surface is defined as Z = 0, and the 
depth considered to be infinite, then 



is the function defining the vertical temperature distribution. The 





( 2 ) 



whose corresponding frequency function is 




dZ Jiff 



(3) 




where CT greater than 0 and 



m are constants, then the frequency function becomes 







(4) 



The function described by (4) has the following characteristics: 



(a) it is symmetric about Z - m 



(b) it has two symmetric points of inflection at 
Z-nat 0* , and 



(c) it has a maximum rate of change 





(5) 



center of the thermocline, or the center of the curve of the 




distribution, is at 
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2f" ! 1 ! I ! I I r ! ^dT 

Normalized Temperature dZ 




Fig.l Gaussian thermocline and its frequency 
function (first derivative) 



k> 



Z =Z C = -m 



( 6 ) 



The top and bottom of the thermocline are then defined as being located 
at the points of the maximum rate of change of slope of theT(3u™j curve 
namely. 



Zt = - nn 1 - cr J3 1 (7) 

z bt - -no - cr JT (8) 



which are illustrated in Figure 1. 

The handling of bathythermograph data leading to the application of 
results (6), (7), and (8) is very simple, since values of m and a may be 
computed directly from the (t,Z) distribution curve. Consider the first 
two moments of the (T,Z) curve about the origin, Z=0, whereby the (nth) 
moment of the curve about the origin is meant the mean values of 
defined as 2.^ > where n is a positive integer. Then with the normal 

curve, the first moment is identical with the center of the thermocline: 



rT(o) 



Zc = 



I 



Z(T) <JT 



(9) 



T(O') -T(-oo) 

where T(0) is the temperature at the surface Z = 0 and T (-oo) is the 
temperature at the bottom of the lower isothermal layer which, in this 
idealized case, is assumed to extend to the bottom of the ocean. The 
second moment 



^T(o) 



7 } - 






Z Z dT 



( 10 ) 



can be used to compute CT , since the variance. O’ , is given by 



<? 





(ID 
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In practice, several simplifications are possible. The temperature T (-go) 
is measured at the bottom of the bathythermograph cast. The differentials, 
dT, become finite differences, ^T, and the integrals are replaced by sum- 
mations. If the depths of the temperatures are read at one degree intervals, 
for example, then the T f s can be replaced by unity. 

One must remember that although integration is a continuous process, 
the less accurate summation over finite intervals is not. The larger the 
interval, the less accurate is the summation. The working value of sigma, 
as computed by the routine, cannot exactly equal the correct value obtained 
from integration over the normal curve. Thus equations (7) and (8) must 
be modified by some weighing factor depending upon the temperature inter- 
val chosen. Allowing this factor to be defined as, sk, equations (7) and 
(8) become 



where a temperature interval of 0.5 degrees gives sk equal to 1.47 and a 
temperature interval equal to 1.0 degrees gives a value for sk equal to 



In summary, the objective definition of the top, bottom and center 
of a Gaussian thermocline requires only readings of (a) surface temperature, 
(b) bottom temperature, and (c) depths at some constant pre-chosen tempera- 
ture interval. 

The center of the Gaussian thermocline is 




(7A) 




(8A) 



1.30. 




( 12 ) 
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Since 




( 13 ) 



Then 



l 



2(TnV - 7J- 



(14) 



The top and bottom of the Gaussian thermocline readily follow from (7A) 
and (8A) . 

One problem remains: the decision as to whether or not a normal 

distribution exists; lacking this distribution, the inaccuracies built 
into this process are not tolerable. This determination can be accomp- 
lished by computation of the third moment. However, in the following 
section it is apparent that this computation is not necessary. 

D. The non-gaussian Thermocline 

The analysis of the more common type of thermal structure, the non- 
gaussian thermocline, requires the computation of finite differences on 
Z(Tn). The first finite difference will be defined as 



Generally, the first finite difference is always positive. The second 
finite difference will be negative until, in this case, the point of 
inflection of the (T,Z) distribution curve is reached, at which point 




(15) 



and the second finite difference as 




(16) 




becomes non-negative. If the first non-negative second 



finite difference is 




then the center of the thermocline is 
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In order to locate Z^, only the upper portion of the (T,Z) distribu- 
tion curve is considered, where a (2.^ Zcl synthetic lower portion to 
this curve segment is constructed so that together with the upper real 
segment a normal distribution is formed* From this curve, 0“ is com- 
puted by methods described in the previous section and is given by 
formula (7A) . (see Figure 2b) . 

Similarly, only the lower (*£ — ~Zc\ part of the (T,Z) distribution 
curve is utilized in computing Z ^ • An upper portion of the curve 

is added to the segment and a new value of & is computed* Z^is 
then calculated by using formula (8A) . (See Figure 2c) . 

E. Irregular Thermal Features . 

1 * Transients 

Transients, according to Boston (1966), are usually considered 
to be less than one degree in magnitude and might not be seen by a greater 
thermal interval. Thus, it is necessary to have a small grid spacing, 
say 0.1 degrees, to identify these transients. Spilhaus (1938) relates 
that the bathythermograph is seldom more accurate than 0.5 degrees, 
borne out by the testing of Stewart (1963), due to such problems as 
dull stylus, hysteresis, and calibration. Even if the instrument were 
exact, and hysteresis were not present, the thickness of the trace limits 
reading accuracy to about 0.2 degrees. Hence, although the grid spacing 
chosen for this program is unrealistic at 0.1 degrees, it is retained with 
the possibility of more accurate thermal instruments, or more accurate 
bathythermograph readings, in the future. 

If such accuracy were possible (see Figure 3) transients could be 
identified by simply identifying the sign changes of . If the 
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a. Complete non-Gauss, 
thermocline 



b. Segments used to 
compute Zt 




c. Segments used to 



compute Zbt 

Fig.2 A non-Gaussian thermocline 
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a. Normal scale 



b. Expanded scale 



Fig. 3 Normal and expanded scale transient 
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first negative value encountered was AZ(Tn') then the beginning of the 
first transient would be at Z(T^-h^ ♦ Then the end of the first thermal 
transient and the beginning of the second transient is signified by a 
sign change of AZ (Tr\) from positive to negative. If this occurs at 

then the end of the first transient occurs at Z(Tn^ and 
the beginning of the next occurs at * This process will con- 

tinue until the top of the first major thermocline is encountered; 
generally we are interested in only the transients in the upper isothermal 
region, 

2 . Multiple the rmoc lines 

The analysis of a thermal distribution containing more than one 
major thermocline proceeds along the lines just discussed in the pre- 
vious section on transients. In Boston’s article (1966) he states that 
a grid spacing of one degree is generally adequate if there are only two 
thermoclines , and a grid spacing of 0.5 is necessary should there be more 
thermoclines present. 

As in the transient analysis the sign changes of AZCTd) are 
examined. If the first positive value is AZCTrO then the center of 

the first thermocline would be located at “Zci • The center of 

the second thermoc line , Zc£ ® i s indicated when the sign of AZXtW) 
changes from negative to positive. If this occurs at &"2-CTV^ , then, 
as in the first thermocline, the center is located at Z(Tnr “ Zc.2. * 

The centers of all subsequent thermoclines are found by continuing this 
process throughout the length of the trace. 

When the centers of the shallowest and deepest thermoclines have been 
identified, the method outlined in Section D, the non-gaussian thermocline 
will be used to determine the top of the shallowest thermocline and the 
bottom of the deepest thermocline. Then considering the fact that c&(jA 
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is a measure of the slope of the (T,Z) curve, the maxima in A Z (Tr^ 
corresponds to the most nearly isothermal conditions. Thus these maxima 
can be used to separate the thermoclines . Once this seperation has been ac- 
complished, the procedure from Section D can be used again to locate the 
top, bottom, and center of each individual thermocline. Considering the 
fact that the thermal structure is assumed to be continuous, with only 
small spacing between consecutive thermoclines, there is little need in 
making the identification any more accurate, or finer, than the subdivi- 
sions that are indicated by the maxima. Figure 4 shows an example of 
real traces in the ocean with this multiple thermocline feature. 

3. Inversions 

Large inversions (+ 0.5 degrees Celcius) are usually found 
where water sources of different origin are close together or where 
strong velocity shears are encountered. Tully and Giovando (1963) docu- 
ment this feature. Uda (1963) and Tully and Dodimead (1963) also show 
that in regions such as the subarctic, where there is excess precipita- 
tion, the summer profile may be characterized by a nearly permanent sum- 
mer inversion located at about 100 meters depth. 

The main feature (see Figure 5) of an inversion is the sign 
change of the slope of the (T,Z) distribution curve. The major points of 
interest to the investigator are usually considered to be: 

(a) upper and lower boundries denoted by Zci and Zc2 
respectively, 

(b) magnitude of the inversion, ZUPP and TUPP, being the 
spacial and thermal magnitudes, respectively, 

(c) minimum temperature, TMIN, 

(d) depth of minimum temperature (depth of axis of sound 
channel), ZMIN. 
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Fig. 4 Triple thermocline - Monterey 
Bay (April, 1968) 
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Fig. 5 Typical BT trace with an Inversion 
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Referring to Figure 5, it should be noted that this system of analysis will 

require a more elaborate notation system, since the curve is no longer 

single valued (that is, more than one depth, Z , will have a specified 

n 

temperature) . In order to distinguish between these depths the notation 
Zj (T n ) will be introduced where the subscript j refers to the order 
in which a multiple occurrence happens from the surface down* Thus, if Tn 
occurs at three depths, the shallowest depth would be indicated asZj(Tn)> 
the intermediate depth Z^On^ , and the deepest as Z^OTV*) . As the 
trace is read vertically from the surface down, the trace appears to be 
single-valued (one depth has one temperature), then double-valued (one 
temperature at three depths), before returning to a single-valued state 
again. Since the finite differences must be computed, its first and 
second finite differences will be defined as &2j(Tn'| and ^ZiCTV,). 
respectively. 

The upper and lower boundaries of the inversion occur at the first 
double-valued temperature. If this temperature is Tn, then Zci — Z^CT/^ 
and Zci = Z 2 Cr n 'i . The last double-valued temperature, or a lower- 
valued single temperature included between the last double-valued tempera- 
tures, is the minimum temperature, TMIN, and the depth of this minimum 
temperature, ZMIN, is defined from • The magnitude of the 

thermal inversion, TUPP, is given by the absolute value of the difference 
between the first and last double-valued temperatures. Once these in- 
version features have been identified, pertinent features above and be- 
low the inversion still must be identified by the methods outlined in 
Section D. 
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3. 



Example of Input 



Two methods of input to this program are data deck and tape deck. 
These two examples illustrate these methods. 

A. Data Deck Input 

In using a data deck input, where an analyzer does not have 
access to the digitizing machines and must key-punch his data deck, it 
is possible to input both positional and curve data. The data deck must 
be made up as follows (see Figure 6) ; 

(1) The data must be printed from the bathythermograph 
curve, at some pre-chosen thermal interval, i.e., - 0.5 degrees. The 
data are punched on the data aperture cards by punching the temperature, 
followed by its depth as picked off at the chosen delta-T. Each data 
bit is allowed ten spaces on the aperture card and no other data can be 
included in this ten space unit, so that a maximum of four data sets 

can be punched per data card,' (since there are only 80 spaces allowed per 
card). The thermal data is continued through as many aperture cards as 
necessary to insure complete recording of the entire trace. 

(2) Immediately following the completion of the first 
trace, the operator must insert M Control Card 1 M , which is a dummy con- 
trol card consisting of "0.0" substituted for the temperature, T, and 
any "non-zero" value substituted for the depth, Z. This dummy card 
signals the main program that the reading routine has completed reading 
the data trace and is ready to read the positional data, if it exists, 
or proceed on with the analysis of the trace. 

(3) Following "Control Card 1" we insert the positional 
data card, if we desire to read in this information. This card contains 
the geographical position, both latitude and longitude to the nearest 
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Trace Data Card 
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Control Card 1 (1st) 
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Fig.6 Card Format and Data Deck Assembly 
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hundredth of a minute, the local date, the Greenwhich Mean Time and date 
of the bathythermograph cast, and the serial number of the bathythermo- 
graph, If the positional data card is used, then every position must be 
filled in exactly as noted in Figure 6 and Appendix II. If the neces- 
sary data is not known, then the unknown spaces must be filled with "0 M, s. 
The format on this card is absolutely vital and inflexible. 

(4) Next we repeat the above procedures (1) through (3) 
for each bathythermograph trace desired. It is not necessary to insert 
a positional card with each trace, but if the positional card is used it 
must be complete, as mentioned above. 

(5) When the entire data package has been made up of trace 
data and positional data cards, then n Control Card 2" is inserted as the 
last card to signal the program that the data is finished and the pro- 
gram can terminate. "Control Card 2" consists of a normal data card for- 
mat, except the "0.0" is inserted in position for both temperature, T, 
and depth, Z. This is always the last card in the data deck (if we do 
not consider the programming control cards used by the system), and will 
always follow the positional data card for the last trace or if this data 
is not used, then it will follow the last "Control Card 1." The entire 
data deck is assembled as shown in Figure 6, and may be of any length. 

The only programming prerequisites are that the thermal inverval be a 
constant for each individual trace, although it is not necessary for all 
the traces to have the same interval, and the input format. Figure 6, 
must be followed rigorously. All other necessary computations, such as 
sigma, weighing factor, etc., are contained or computed within the 
program. 
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B. Magnetic Tape Data 

When the input data can be digitized, by such facilities as the Fleet 
Numerical Weather Central, Monterey, or the Canadian Oceanographic Data 
Center, a few modifications must be made to the basic program reading meth- 
ods as they presently exist in Appendix II* In order to use magnetic tape 
input it is necessary to remove the cards that are serialized (right hand 
margin) 0185 - 0204 and insert in their place new cards designed to read 
the magnetic tape* It is important to note here that the format of the 
tapes and even the type of tapes makes a vast difference in the make up 
of the reading cards and great care must be used in inserting the cor- 
rect control and reading cards in order to insure not shutting down the 
computer. Before using tapes, be sure to check the individual instruc- 
tion manual and the individual tape specifications and rigorously follow 
the instructions on the format for the make up of these new cards to be 
inserted. The following information must be considered in each case: 

(a) whether seven or nine track tape, 

(b) data density, 

(c) whether even or odd parity, 

(d) whether the data are translated or not, 

(e) if translated, what machine language, 

(f) record length, 

(g) whether or not the data are blocked, 

(h) if blocked, what size* 

All this information is vital in making up the correct system control 
cards for the tape reading unit. Once the system control cards are cor- 
rectly made up and inserted around the program, then a normal ’'read" card, 
for example; 
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READ (8,1, END=60) 

1 FORMAT(8A5) 

is used, where, in this example, we are reading on tape reader "8", 
sixty traces consisting of eight data sets per trace. 
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4 . Computer Program Operation , 



The computer first reads in the data for an individual trace, as 
outlined in Section 3, and then reads in the positional data (if existent 
in data deck) * Once the individual trace is read into the program, the 
computer takes over, using the analysis method programmed by this routine 
to start the investigations for the features outlined in Section 2„ Be- 
fore starting any further investigations the program looks at the thermal 
interval and choses one of three possible values for the weighing factor, 
SK. This value is then carried throughout the calculations for this trace,, 
If the positional data is included, the first output is a formalized print- 
out of this informational heading. Once this has been executed, the main 
program starts a searching routine to identify various possible features 
before calling the specific sub-routines necessary to analyze this data. 

The program first looks for double-valued temperatures in the vertical 
dimension. If these double-values exist the program calls the inversion 
subroutine, INVERS, to calculate the top and bottom of the inversion and 
its spacial and thermal magnitude. This subroutine first calls a special- 
ized subroutine, 0RDER1 , that takes all the input data from that trace 
and checks to insure that the data is in correct order from the surface 
down (consecutive readings) and if it is not, reorders the data in the 
correct sequence. Once the data is in correct order, then the program 
does a quick safety check to insure that it does not have a non- inversion, 
(where the temperatures first increases with depth and then decreases, 
giving not an inversion, but a positive surface gradient followed by the 
normal main thermocline, a feature that the main program might identify 
as an inversion). If this condition does exist, the subroutine terminates 
and returns the flow to the main program for further investigations. If 
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an inversion exists ,then, by the use of do-loops, INVERS calculates the 
top, bottom, and magnitudes of the inversion as outlined in Section E.3. 
Once the data has been analyzed and these results are available, the 
subroutine then calls four specialized printing subroutines, PRINT3 
through PRINT6, to print out the data in a useful format. INVERS then 
returns the order flow to the MAIN program for further investigations. 

The MAIN program then investigates the thermal interval, delta-T, 
to insure that if the interval is less than 0.2 degrees a thorough 
investigation for the presence of transients is accomplished. If the 
interval is 0.2 or less, the MAIN routine calls the transient investiga- 
tion routine, TRANS, which proceeds to analyze for the presence of transi- 
ents as outlined in E.l. Upon completing the analysis, the results are 
printed out, if any exist, and the flow returns to the MAIN routine. 

Once returned to the MAIN program, the final investigation is made 
into the possibility of multiple thermoclines . The MAIN program calls 
the multiple thermocline subroutine, MANYT, which utilizes analysis 
procedures outlined in Section E.2. Within this subroutine is the pro- 
vision to go to another subroutine, NONGA, designed to handle the more 
specific case, should only a single main thermocline exist. The NONGA 
subroutine uses the procedure outlined in Section D to find the top, 
center, and bottom of the thermocline. In both subroutines, looking at 
thermoclines, the results are printed out by utilizing the specialized 
printing subroutines, PRINT1 and PRINT2. 

A flow chart (simplified) for the entire process is shown in Ap- 
pendix I, and the entire program can be inspected by referring to Ap- 
pendix II. It is worthy of note, that this program will not identify all 
three features, transients, multiple (or single) thermoclines, and in- 
versions, in a single trace as the existance of all three would be very 
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rare. However, the program will look for transients in the upper near- 
isothermal layer and then investigate the possibility of an inversion be- 
low this area, or if the inversion does not exist then look at the thermo- 
cline(s) in the region below the near isothermal layer. 

Also, as an additional feature the program can be called to print 
out a graphic display of the input data (Appendix III), thus giving a 
visual representation of the data that is printed out on the output data 
sheets. This auxiliary feature can be called by inserting the following 
cards in the MAIN program immediately following the card serialized 0199 . 

CALL SEE 13 (T 5 Z s ICOUNT) 

which will call the SEE13 subroutine that utilizes the library pro- 
gram, "DRAW", to give graphic display of the input data set and fits a 
curve to these points. Before using this call, the user must exercise 
care and carefully read the rules for usage of "DRAW" in subroutine 
n SEE13 M as are outlined in Appendix III. 

Another feature included as part of this thesis package, is a pro- 
visional set of cards (Appendix IV) that will convert the program as 
contained in Appendix II into the slightly different machine language 
required when the program is used on the CDC 6500 model computers, in- 
stead of the IBM 630 models for which this program is written. Before 
using this conversion program two steps must be taken, 

(1) all data and system control cards must be removed 
and the remaining portions of the included program 
become the data deck for "BCD 90 

(2) All output formats that involve printed words or 
combinations of printed words and numbers must be 
converted to "H" format. 
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5. 



Output from the Computer. 



The output from the computer consists of several possible combina- 
tions of data.,. Included are examples of a few of the possible combina- 

tions. All the output from a given trace is printed on a single indivi- 
dual page, headed by the trace number within the data sequence and the 
positional data for that slide (if any). Figure 7 shows one of the sets 
of real data that was taken from some data traces made by Tully in the 
vicinity of Ocean Station "Papa” in 1963. Figure 7 is a hand-drawn en- 
largement of a real test trace containing both transients and a single 
main thermocline, typical of that region, used in testing out this pro- 
gram. Figure 8 is a typical output from the program's analysis of this 
trace and lists the positional data (which was partially available) and 
the date and time of the trace. The program then reveals the existence 
of two transients, by printing out the location of the start and end of 
each transient. It also gives the information on the single major thermo- 
cline. In this case the data was read in from a data dekc after being 

transcribed from the trace onto aperture cards. The data interval was 
taken in this case to be 0.1 degrees. Figure 9 shows the typical print- 
out of a hand-constructed case, also used to test program, centered 
around a thermal inversion (Boston 1966). Figure 10 is the data output 
for this created case. Figure 12 is an example of the printed output 
that would be expected from a trace (see Figure 11) containing multiple 
thermoclines . In this case two major thermoclines exist in the trace. 

For interest. Figure 13 shows the output from the drawing routine called 
from SEE13. It can be seen that computer drawn trace of Figure 13 is an 
adequate enough representation of Figure 12 to enable the user to use 
this type of computer output as a rough aid in insuring correct working 
of the computer program or other rough calculations where a visual presenta- 
tion might be desirable. 
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Fig. 7 Ocean Station 'Papa' — BTTrace 

23 April 1963 
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START OF INVESTIGATIONS IN BATHYTHERMOGRAPH TRACE NUMBER 4 

GEOGRAPHIC POSITION OF BATHYTHERMOGRAPH TRACE: 

LATITUDE 36-36.1 NORTH 

LONGITUDE 113-44.7 WEST 

DATE AND TIME OF BATHYTHERMOGRAPH TRACE: 

DATE .10/23/63 

TIME .100832 

BATHYTHERMOGRAPH SERIAL NUMBER: 

BT SERIAL NO 663977 

NO TEMPERATURE INVERSIONS EXIST IN THIS TRACE. 

TRANSIENT DATA AS COMPUTED BY PROGRAM: 

BEGINNING OF FIRST TRANSIENT 21.0 

END OF TRANSIENT 55.0 

BEGINNING OF NEXT TRANSIENT 60.0 

END OF TRANSIENT 87.2 

MULTIPLE THERMOCLINE DATA AS COMPUTED BY PROGRAM: 

CENTER NUMBER 1 LOCATED AT . 148.0 

CENTER OF THERMOCLINE 149.6 

TOP OF THERMOCLINE 99.3 

BOTTOM OF THERMOCLINE 196.8 

END OF DATA FOR THIS BATHYTHERMOGRAPH SLIDE. 

NUMBER OF DATA PAIRS ANALYZED - 18 

Fig. 8 Example of Output for Real Data Curve 
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Fig. 9 Typical Inversion 
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START OF INVESTIGATIONS IN BATHYTHEKMOGRAM TRACE NUMBER 27 

GEOGRAPHIC POSITION OF BATHYTHERMOGRAPH TRACE: 

LATITUDE 00.00 

LONGITUDE 00.00 

DATE AND TIME OF BATHYTHERMOGRAPH TRACE: 

DATE 00/00/00 

TIME 000000 

BATHYTHERMOGRAPH SERIAL NUMBER: 

BT SERIAL NO 000000 

THE FOLLOWING INDICATES THE COMPUTED TOP OF THE TEMPERATURE 
INVERSION, BY NON-GAUSSIAN THEODS: 

TEMPERATURE INVERSION TOP 66.0 

THE FOLLOWING INDICATES THE COMPUTED BASE OF THE TEMPERATURE 
INVERSION, BY NON-GAUSSIAN METHODS: 

TEMPERATURE INVERSION BASE 160.1 

THE FOLLOWING FIGURES DEFINE THE MINIMUM TEMPERATURE OF THE 
TEMPERATURE INVERSION: 

MINIMUM TEMPERATURE 11.0 

DEPTH OF MINIMUM TEMPERATURE 139.4 

THE FOLLOWING FIGURES DENOTE THE TEMPERATURE AND SPACIAL 
MAGNITUDE OF THE TEMPERATURE INVERSION: 

THERMAL MAGNITUDE 4.0 

SPACIAL MAGNITUDE 94.1 

NO THERMAL TRANSIENTS EXIST IN THIS TRACE. 

NO MULTIPLE THERMOCLINES EXIST IN THIS TRACE. 

END OF DATA FOR THIS BATHYTHERMOGRAPH SLIDE. 

NUMBER OF DATA PAIRS ANALYZED - 

Fig. 10 Example of Output for Inversion 
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Fig.11 Hand drawn Bathythemogram 
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START OF INVESTIGATIONS IN BATHYTHERMOGRAM TRACE NUMBER 17 

GEOGRAPHIC POSITION OF BATHYTHERMOGRAPH TRACE: 

LATITUDE 000.00 NORTH 

LONGITUDE 000.00 WEST 

DATE AND TIME OF BATHYTHERMOGRAPH TRACE: 

DATE 00/00/00 

TIME 00000 

BATHYTHERMOGRAPH SERIAL NUMBER: 

BT SERIAL NO 000000 

NO TEMPERATURE INVERSION EXIST IN THIS TRACE. 

NO THERMAL TRANSIENTS EXIST IN THIS TRACE. 

MULTIPLE THERMOVLINE DATA AS COMPUTED BY PROGRAM: 

CENTER NUMBER 1 LOCATED AT ......... 72.0 

CENTER NUMBER 2 LOCATED AT 144.0 

CENTER OF GAUSSIAN THERMOCLINE 75.0 

TOP OF GAUSSIAN THERMOCLINE 53.5 

BOTTOM OF GAUSSIAN THERMOCLINE 94.4 

SEPERATION #1 BETWEEN THERMOCLINES 100.0 

CENTER OF GAUSSIAN THERMOCLINE 138.0 

TOP OF GAUSSIAN THERMOCLINE 104.9 

BOTTOM OF GAUSSIAN THERMOCLINE ..... 152.6 

END OF DATA FOR THIS BATHYTHERMOGRAPH SLIDE. 

NUMBER OF DATA PAIRS ANALYZED - 8 

Fig. 12 Example of Output for Multiple Thermoclines 
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TEMPERATURE AGAINST DEPTH (Z POS. UP) 

FIG. 13 
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All the preceeding traces were read in in data deck format, with 
varied intervals of delta-T, 0.1 degrees in the first example, 1.0 
degrees in the second case, and 0.5 degrees in the latter case. There 
are other possible combinations, but these three cases outline the general 
format and the answers that can be expected. The other cases are just 
various combinations or omissions of the above three cases. 
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6 . 



Discussion and Conclusions Concerning Results . 



Although no error analysis was undertaken for this computer pro- 
cess, it is obvious that there are some errors built into the program 
that could perhaps be reduced by further Investigation and refinement. 

In several segments we use a Gaussian based method and thus there is 
some error introduced in the computation of the top and bottom of the 
thermocline due to having to compute a sigma in the program. Also, in 
this program the multiplication factor, SK, is calculated only roughly 
based on the thermal interval. This factor can be improved, but cannot 
be exact in the computer, due to truncation error; besides, much storage 
space is needed to store the large table of values applicable where each 
tabular figure is only good for one specific set of data. 

There are certain errors inherent in the computer, but since we do 
not require great precision, these errors can be ignored due to the fact 
that the bathythermograph itself is not a precise instrument. These 
errors include truncation and round-off errors in every computation. 

Since we are requiring printout to only one decimal place it is doubt- 
ful that these errors will ever be hoticable. Also there are errors 
built into the functions that are interral to the computer, such as the 
"square root," "multiplication," and "addition" routines, but again these 
will not be significant in the accuracy that we require. 

There are certain advantages to the use of data tapes versus magne- 
tic cards, when great volumes of data are to be processed. The time 
required to punch the data cards and to check the results is a function 
of the experience of the operator. The computer time is approximately 
25 seconds for card input computations, plus additional time if the draw- 
ing of graphs is required (approximately three minutes per graph). This 
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time can be cut approximately in half if the program compiled by the 
computer is put on tape and data input is in that form, as the compiler 
takes approximately half the time now used in making a computer run. 

There are other equipment errors, that as yet have not been evalu- 
ated, that are concerned with the input data, and are still vital to the 
person desiring to analyze the traces by use of this program. Such items 
as hysteresis, where a double trace is formed on the bathythermograph 
slide, due to slight differences in the temperature when raising the 
Bathythermograph as compared to when lowering it cause inaccuracies* 
dull styluses can give traces that are too wide to be useful. Damaged 
equipment that gives unacceptable temperature fluctuations, such as 
variations from bucket thermometer temperatures, or other signs that the 
bathythermogram is inaccurate although the trace itself appears accept- 
able is another input error. Also there are errors inherent in the digiti- 
zer that can vary depending upon the equipment used and the usual problems 
inherent to electronic curve followers or other devices. 

The results of this program, in the fifty test cases run, appear to 
give slight errors as compared to the hand analysis of the same data. 

Ten of the test cases were hand constructed curves, designed to test 
the program for specific features or combinations of these features, and 
the remaining data were real data obtained from the Fleet Numerical Weather 
Central or taken from other written works, such as Tully and Dodimead 
(1957) and Pattullo and Cochrane (1951), and a few traces taken in 
Monterey Bay. Both tape inputs from Fleet Numerical Weather Central 
and punched tape data decks were used. There were no apparent differ- 
ences between the outputs from the digitized data and the data on 
aperature cards, as might be expected assuming that care was taken by 
the person picking off data points from the curves. The largest error 
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in the tested cases, was found in a very steep thermocline, real data 
case, where the depth error was 4.5 feet for the top of the thermocline 
and 3.0 feet for the base of the thermocline in a thermocline that was 76 
feet in vertical magnitude. In this case the center of the thermocline 
was only 0.3 feet from the hand calculated value. The thermal interval 
was 1.0 degree Farenheit and the thermocline in question was analyzed 
by the Gaussian technique. 

Thus it can be said that this program accomplishes, with acceptable 
accuracy, considering the analytic methods used, in seconds what a well 
trained operator would take tens of minutes to accomplish. 

This program can be expanded, by the inclusion of additional sub- 
routines, and their respective "call" cards in the MAIN routine without 
greatly increasing the time needed in computing the desired results. 
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7. 



Recommendations . 



As previously stated, there are other features that can be added 
to this program by merely including additional subroutines. Additionally, 
this same programming technique can be used, with a few basic changes, to 
analyze other oceanographic data curves. Before accomplishing any con- 
versions of this program to handle other data, the user must first look 
at the basic form of the traces he plans to analyze. Traces, other than 
temperature, may not be describable using the same parameters as des- 
cribed by Tully (1964). Such items as the magnitude of the fluctuations 
of other quantities must be examined as they may exceed the number that 
this program can handle. For example, it may be common, in some other 
oceanographic parameter, to have more than three stable features corres- 
ponding to the thermoclines in this program, and thus the program would 
have to be modified to allow for additional curve computations. Another 
correction that would have to be made would be the rewording of the print- 
out statements. 

Further investigations should be undertaken, to include the improve- 
ment in accuracy of sigma and SK parameters and the possibility that such 
curves as the hyperbolic tangent or Chi-squared curves might give a better 
curve fit. 

As this program has yet to be run on the CDC 6500 computer of the 
Fleet Numerical Weather Central, this should be undertaken as soon as 
possible. For this type of program, the speed of computation on the CDC 
6500 computer is approximately three times that of the IBM 360. This 
CDC 6500 computer also contains more storage space than the IBM 360 
computer and thus longer tapes or card inputs may be used. An additional 
advantage of this increased storage is that investigation could be shifted 
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from the present two dimensional set up (temperature versus depth) to 
a three dimensional investigation (multiple traces analyzed simultaneously) 
giving a field representation instead of just specific points. 
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PRINT OUT 

Print out; '"NORTH* latitude and "EAST" 
longitude and also date, time, and serial 
of bathythermograph. 







PRINT OUT 






Print out; "NORTH" latitude and 'VEST" 
longitude and also date, time, and serial 
of bathythermograph. 












PRINT OUT 

Print out; "SOUTH" latitude and "EAST" 
longitude and also date, time, and serial 
of bathythermograph. 
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Subroutine TRANS 



Inputs from: MAIN 

Z(I), T(I), ICOUNT 



Set do-loop counters, based on size of 
data set, ICOUNT. 



Bydo-loop compute first difference, 
Z6(N) , on depth. 



By do-loop, compute second difference, 
Z 7 (141) , on Z6(N) . 



-YES _ 



By loop, check sign of second difference 
to see if it is greater than zero. 



PRINT OUT 

Print out, "Beginning of first transient" 



YES 



NN.LE.M 



CONTINUE 



NN.LE.M 




YES 
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Set up localization loop to find extent of 
the lower portion of the thermocline 
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PRINT OUT 



Print out "Base of thermocline located at 
depth . " 
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APPENDIX II 

OBJECTIVE THERMOCLINE ANALYSIS PROGRAM 
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WRITE-UP OF "DRAW" SUBROUTINE 
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LABFl IF 4 CURVE IS BRING DRAWN <ITYPF=0>, t_AB c L IS A <=01)* 

CHARACTER BCD CURVE IDRNTIPIFR WHICH WILL BE RRRRODUCF 0 
AT THE END D F THR CURVE. RPR EXAMPLE » LABEL EQUALS 



AH INF. I c POINTS Aoe RFIND PLOTTED, L A R p I IS AN R- 
CHARACTER I DFN T I F I C R . THF FIRST FOtJR CHARACTERS AR C 
REPRODUCED WITH the ejRST omTTFD POINT, AND TH C L A c T 
FOUR WITH thf LAST POINT. S C T TO RIAN‘S ANY UNWANTED 
CHARACTERS. 
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"BCD" CONVERSION PROGRAM 
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